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A detection or nondetection of primordial non-Gaussianity in the CMB data is es- 
sential not only to test alternative models of the physics of the early universe but also to 
discriminate among classes of inflationary models. Given this far reaching consequences 
of such a non-Gaussianity detection for our understanding of the physics of the early uni- 
verse, it is important to employ alternative indicators in order to have further information 
about the Gaussianity features of CMB that may be helpful for identifying their origins. 
In this way, a considerable effort has recently gone into the design of non-Gaussianity 
indicators, and in their application in the search for deviation from Gaussianity in the 
CMB data. Recently we have proposed two new large-angle non-Gaussianity indicators 
which provide measures of the departure from Gaussianity on large angular scales. We 
have used these indicators to carry out analyses of Gaussianity of the single frequency 
bands and of the available foreground-reduced five-year maps with and without the KQ75 
mask. Here we extend and complement these studies by performing a new analysis of 
deviation from Gaussianity of the three-year harmonic ILC (HILC) foreground-reduced 
full-sky and KQ75 masked maps obtained from WMAP data. We show that this full- 
sky foreground-reduced maps presents a significant deviation from Gaussianity, which 
is brought down to a level of consistency with Gaussianity when the KQ75 mask is 
employed. 

Keywords: Gaussianity; cosmic microwave background, inflation, physics of the early 
universe. 



1. Introduction 

A detection or nondetection of primordial non-Gaussianity in the CMB data is cru- 
cial not only to discriminate inflationary models but also to test some alternative 
scenarios for the physics of the early universe. However, the extraction of primor- 
dial non-Gaussianity is a difhcult enterprise since several effects of non-primordial 
nature can produce non-Gaussianity in the CMB data. Clearly the study of de- 
tectable non-Gaussianities in the WMAP data must take into account that they 
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may have non-cosmological origins as, for example, unsubtracted fore groun d con- 
tamination, unconsidered point sources emission and systematic errors P^ISI Devia- 
tion from Gaussianity may also have a cosmic topology origin (see, e.g., the review 
articles Refs. |4]and related Refs. [5]). If, on the one hand, different statistical tools 
can in principle provide information about distinct forms of non-Gaussianity, on the 
other hand one does not expect that a single statistical estimator can be sensitive 
to all possible forms of non-Gaussianity in CMB data. It is therefore important 
to test CMB data for Gaussianity by using different statistical indicators to shed 
some light on its possible causes. In view of this, a great deal of effort has recently 
gone into verifying the existence of non-Gaussianity by employing several statistical 
estimators.!^ 

Recently have we proposeci^ two new large-angle non-Gaussianity indicators, 
based on skewness and kurtosis of large-angle patches of CMB maps, which provide 
measures of the departure from Gaussianity on large angular scales. We used these 
indicators to search for the large-angle deviation from Gaussianity in the three and 
five-year single frequency K, Ka, Q, V, and W maps with and without a KQ75 
mask. We have found strong deviation from Gaussianity in the unmasked maps, 
whereas a KQ75 mask lowers significantly the level of non-Gaussianity (for details 
see Ref. [7|. 

Motivated by the fact that most of Gaussianity analyses with Wilkinson Mi- 
crowave Anisotropy Probe (WMAP) data have been carried out by using CMB 
frequency bands masked maps, and that sky cut can in principle induce bias in Gaus- 
sianity analyses, in a more recent papeJ^ we have carried out an analys is of Ga us- 
sianity of the available full-sky foreground-reduced five-year CMB maps l ^ l ^^ l ^^ l by 
using the statistical indicators of Ref. [3 

We have shown that the available full-sky five-year foreground-reduced maps 
present a significant deviation from Gaussianity, which varies with the foreground- 
cleaning procedures. We have also shown that there is a substantial reduction in 
the level of deviation from Gaussianity in these full sky maps when a KQ75 mask 
is used. Our main aim here is to extend and complement our previous worlPl by 
performing a similar analysis of Gaussianity of the three-year harmonic ILC (HILC) 
mapJi^ foreground-reduced full-sky and KQ75 masked maps. To this end, in the 
next section we give an account of the large-angle non-Gaussianity indicators of 
Ref. 7, while in the last section we apply our indicators to perform a Gaussianity 
analysis of the HILC three-year full-sky and KQ75 cut-sky maps, and present our 
main results. Our principal conclusion is that the HILC full-sky foreground-reduced 
maps presents a significant deviation from Gaussianity, which is reduced to a level 
of consistency with Gaussianity when the KQ75 mask is employed. 

2. Non-Gaussianity Indicators and Associated Maps 

A constructive way of defining our non-Gaussianity indicators S and K (discrete 
functions defined on S*^) from CMB data can be formalized through the following 
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i. Take a finite set of points {j = l,...,A^c} homogeneously distributed on the 
CMB celestial sphere as the centers of spherical caps of a given aperture 7; 
and calculate for each cap j the skewness and kurtosis given, respectively, by 



where A'p is the number of pixels in the j cap, Ti is the temperature at the i 
pixel, Tj is the CMB mean temperature of the j cap, and a is the standard 
deviation. Clearly, the numbers Sj and Kj obtained in this way for each cap can 
be viewed as a measure of non-Gaussianity in the direction of the center of the 
cap {Oj,(j)j). 

iii. Patching together the Sj and Kj values for each spherical cap, one obtains our 
indicators, i.e., discrete functions S = S{9, (j)) and K = K{9, (j)) defined over the 
celestial sphere, which can be used to measure the deviation from Gaussianity 
as a function of the angular coordinates {0,(j>)- The Mollweide projection of 
skewness and kurtosis functions S = S{d, (j)) and K = K{9, (p) are nothing but 
skewness and kurtosis maps, hereafter referred to them as S'— map and X— map, 
respectively. 

Clearly, the discrete functions S = S{d, (j)) and K = K{9, (p) can be expanded 
into their spherical harmonics in order to determine their power spectra Sg and Kg. 
Thus, for example, for the skewness one has S{9, (f>) — X^fco '^m=-i Yim[9, (f>) 
and Si — {2£ + 1)"-^ l^^mP- Similar expressions obviously hold for the kurtosis 



3. Main Results and Conclusions 

In this section we shall report the results of our Gaussianity analysis performed with 
S = S{9,(j)) and K = K{9,(j)) indicators calculated from the foreground reduced 
HILC full-sky and KQ75 masked maps computed from three- year WMAP data0 
To minimize the statistical noise, in the calculations of S*— map and K—nisq) from 
the HILC foreground-reduced three-year map, we have scanned the celestial sphere 
with 12 288 spherical caps of aperture 7 = 90°, centered at points homogeneously 
generated on the two-sphere by using 

HEALPijPl. 

Figure [1] shows the differential power spectrum of the skewness Sg (left panel) 
and kurtosis Kg (right panel) indicators for i = 1, • • ■ , 10 , calculated from full-sky 



^For a detailed discussion of the indicator briefly presented here we refer the readers to Ref. [7]and 
Ref.ll] 

''We note that in the analysis of Gaussianity with the KQ75 masked maps the implementation of 
the mask is made by removing the pixels inside the masked regions from the set of pixels of the 
each cap whose intersection with the mask is not empty. Thus, the values Sj and Kj for a j*^ cap 
(with pixels in the mask region) are calculated with small number Np of pixels. 




(1) 



K = K{9,cP). 
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S. HILC-3yr: full-sky & with KQ75 



A - HILC-3yr full-sky 
• = HILC-3yr + KQ75 



1 25456789 10 
multipoles (0 



K, HILC-3yr: full-sky & with KQ75 



A = HILC-3yr full-sky 
• = HILC-3yr + K075 



1 25456789 10 

multipoles (0 



Fig. 1. Low £ differential power spectra of skewness \Si — Si\ (left) and kurtosis (right) \Ki — 
Ki\ calculated from the foreground-reduced HILC full-sky and KQ75 masked maps. The 95% 
confidence level (obtained from Monte-Carlo Gaussian maps) is indicated by the dashed line. 

and KQ75 cut-sky three- year foreground-reduced HILC maps. The 95% confidence 
level, obtained from S and K maps calculated from Monte-Carlo (MC) statistically 
Gaussian CMB maps, is indicated in this figure@ To the extent that the deviations 
\S£ — Se\ and \Ki ~ Ki\ for these maps are not within 95% of the mean MC value, 
Fig.[T]shows an important deviation from Gaussianity in full-sky foreground-reduced 
HILC three-year map. This figure also shows a significant reduction in the level of 
large-angle deviation from Gaussianity when the KQ75 mask is used. 

To have an overall assessment power spectra and Ki calculated from the HILC 
foreground-reduced three- year full and cut map, we have performed a test to find 
out the goodness of fit for Sg and Kg multipole values as compared to the expected 
multipole values obtained from S and K maps calculated from Monte-Carlo (MC) 
statistically Gaussian CMB maps. This gives a number for each case that quantifies 
collectively the deviation from Gaussianity. For the power spectra Se and Ki we 
found the values given in Table [T] for the ratio x^/dof (dof stands for degrees of 
freedom) for the power spectra calculated from three-year HILC foreground-reduced 
full-sky and cut-sky maps. 



Table 1. x test goodness of fit for Si and Ki 
calculated from the HILC full-sky and cut-sky 
three- year maps as compared with the expected 
values Si and Ki obtained from MC maps . 



Map 



X-' for Se X for Ke 



HILC full-sky 
HILC KQ75 cut-sky 



1.6 X 10^ 8.8 X 10^ 
0.8 1.5 



Clearly, the greater is the values for x^/dof the smaller are the probabilities, 
that is the probability that the power spectra Si and Ki and the expected MC 
power spectra agree. Thus, from Table [T] is one concludes that the HILC presents 



"^For details on the calculation of these (data and MC) maps and the associated power spectra we 
refer the readers to Ref . UI and Ref. |8] 
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the substantial level of deviation from Gaussianity as detected by the indicators, 
which is reduced to a level that can be considered consistent with Gaussianity when 
the KQ75 mask is employed. 

Finally we note that the relative deviation of the full-sky power spectrum from 
the cut-sky spectrum can be calculated with no reference to the Gaussian MC 
spectra. To this end, we have performed a test to find out the goodness of fit for 
Si and Kg multipole values for the full-sky maps as compared to the corresponding 
cut-sky values. For this relative assessment of power spectra and Kg we have 
found that x^/dof are 1.4 x 10^ and 7.2 x 10^. These values make apparent the 
significant effect of the mask in the reduction of the deviation from Gaussianity in 
the full-sky HILC three- year map, and give information on reliability of the HILC 
full-sky foreground-reduced three-year map as Gaussian reconstruction of the whole 
CMB sky. 
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